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Recent observations of the deeply embedded L 1551 IRS 5 system permit the 
detailed examination of the properties of both the stellar binary and the binary 
jet. For the individual components of the stellar binary, we determine their 
masses, mass accretion rates, effective temperatures and luminosities. For the 
atomic wind/jet flow, we determine the mass loss rate, yielding observationally 
determined values of the ratio of the mass loss to the mass accretion rate, /. For 
the X-ray emitting region in the northern jet, we have obtained the jet-velocity 
and derive the extinction and the densities on different spatial scales. Examining 
the observational evidence within the framework of the x-wind theory leads us 
to conclude that these models are indeed potentially able to account for the 
observational data for this deeply embedded source. 

Subject headings: ISM: jets and outflows — ISM: individual (L1551 IRS 5 jets) 
— stars: formation — stars: pre-main-sequence — binaries: general — accretion, 
accretion disks 
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1. Introduction 

Since their discovery nearly three decades ago, the unexpected phenomenon of outflows 
in star forming regions has remained essentially unexplained. In particular, the processes 
responsible for the acceleration and collimation of the flows present one of the major un- 
solved problems of modern astrophysics (Lada 1985; Shu et al. 1987, 2000; Eisloffel et al. 
2000; Konigl & Pudritz 2000). Considerable progress has been made in the theoretical field, 
but observational results had generally been of too poor quality to make direct comparisons 
meaningful. In this paper, we examine recent various observational results for the deeply 
embedded source L 1551 IRS 5, which, when put together, finally permit the detailed com- 
parison with theoretical models. Optically visible young stellar objects (lower mass loss and 
accretion rates) have previously been addressed by, e.g., Shang et al. (2002, and references 
therein). 

Over the years, IRS 5 has enjoyed a variable status of 'the archetypical CO outflow 
source' and that of 'a pathological case' (e.g., Padman et al. 1997). More recently, renewed 
interest in this object has arisen, in part due to the recognition of its duplicity, as bina- 
rity seems to be a very frequent phenomenon among powerful outflow drivers (M. Barsony, 
private communication). If protostellar binarity is indeed intimately related to the physics 
of generating bipolar outflows, then one obviously wishes to understand the physics of the 
binary itself. We devoted therefore quite some effort to arrive at the understanding of the 
IRS 5 system. Based on the theory of protostellar structure and evolution, we find a solu- 
tion capable of explaining observed characteristics of both IRS 5 itself and its associated jet 
flows, and which, when applied to x-wind models, successfully recovers the physical proper- 
ties of the jets. Thus, in this particular case, a direct dependence on source binarity is not 
evidenced, as the theory of x-winds has been developed for single stellar objects. 

The organisation of this paper is as follows: In Sect. 2, we review the evidence based on 
recent observations in both the optical, infrared and radio spectral regimes. In Sect. 3, we 
use published models of protostellar structure and evolution to derive the physical properties 
of the individual binary components and apply x-wind theory to derive some parameters 
relevant in the present context. In Sect. 4, we discuss our results and, finally, in Sect. 5, we 
briefly summarise our main conclusions. 
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2. The Observational Evidence 

2.1. Mass loss rates from the large scale atomic flow 

Hollenbach (1985) proposed that the luminosity of the [Ol]63/im line can be used 
to estimate the wind mass loss rate. At the distance of 150 pc and for the observed line 
luminosity with an 86" beam 1 (White et al. 2000), this method would yield a mass loss rate 
of 8 x 10" 7 M© yr" 1 . In contrast, based on H I 21 cm observations of the redshifted gas with 
the VLA (Very Large Array, ~ 57" beam), Giovanardi et al. (2000) determined the mass loss 
rate from IRS 5 for the atomic wind as M w>re( j = 1 x lO _6 M yr _1 . Given the momentum 
symmetry of the red- and blueshifted flows (see their Table 2), a mass loss rate of about 
2 x 10~ 6 M yr" 1 seems indicated, a factor of three higher than that based on the [O i] 63 /im 
line. 

The reason for this discrepancy is not entirely clear. One possibility could be that the 
condition, in which the Hollenbach relation is thought to hold, is not met by the flow from 
IRS 5 ( viz. $h ^ TT-o, 5 t's,ioo; where $h is the particle flux in cm s through the shock front, 
and when the pre-shock density n is expressed in units of 10 5 cm -3 and the shock velocity 
v s is in units of 100 km s" 1 ). To some degree this might be supported by the fact that IRS 5 
falls out completely, by more than one order of magnitude, of the empirical I/forbidden iinc — £boi 
relation found for a large number of jet-sources (Edwards et al. 1993; Liseau et al. 1997). 
However, most of these sources, being of T Tauri type, are presumably at a later stage of 
development. 

Source variability could be another, perhaps more likely, reason and in which case 
neither of the two M w estimates would be entirely correct, since both methods are based on 
the assumption of stationarity. Intensity variations from the shocked gas on small spatial 
and temporal scales (< 10", < 0.5 yr) have frequently been observed (Cameron & Liseau 
1990; Liseau et al. 1996) and the [O i] 63 /mi cooling time is likely of this order. As discussed 
by Giovanardi et al. (2000), the HI data pertain to time scales of the order of 30 yr and the 
memory of small scale variations might become 'ironed out' in the global flow. 



1 This refers to the flux measured by the ISO-LWS (Infrared Space Observatory Long Wavelength Spec- 
trometer) toward IRS 5. In addition, another ten positions were observed toward the two CO outflow lobes. 
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2.2. The small scale flows 

2.2.1. The velocities of the optical jets 

On the sub-arcsec to ten arcsec scale, two independent jets from IRS 5 have been iden- 
tified in the optical and in the near infrared (Fridlund & Liseau 1998; Itoh et al. 2000; 
Fridlund et al. 2004). These have been designated as the 'northern' and the 'southern' 
jet, respectively. The two jets distinguish themselves by very different emissivity and ve- 
locity characteristics. The intensity of the northen jet is much stronger and the material 
exhibits much larger velocities. At the projected jet position 2" (300 AU) distant from IRS 5, 
Fridlund et al. (2004) have recently obtained a maximum (blue-shifted) radial velocity of 
maxi) ra j = vcos8 = 430kms _1 . This flow velocity has been maintained at this high level 
over the years (cf. Stocke et al. 1988). Radial velocity and proper motion data imply that 
the jet is moving at non-zero angles with respect to the line of sight, i.e. 6*n > (Liseau 
& Sandell 1986; Stocke et al. 1988; Fridlund & Liseau 1994; Lucas & Roche 1996), so that 
the jet velocity is strictly higher than 430kms _1 , consistent with the shock generated X-ray 
luminosity observed from this jet position (Favata et al. 2002; Bally et al. 2003; Fridlund et 
al. 2004). In contrast, the radial velocity observed for the southern jet is at most 65kms _1 
(Fridlund & Liseau 1998; Hartigan et al. 2000; Pyo et al. 2002; Fridlund et al. 2004) and 
any non-zero proper motion is below detectability in the multi-epoch data of Fridlund et al. 
(2004). 

2.2.2. The densities of the northern optical jet 

At the position of observed highest radial velocity, Fridlund et al. (2004) have estimated 
local volume densities from spectrophotometric observations of lines of [SII], viz. n e = 
(1.8 ± 0.4) x 10 3 cm -3 . This value represents the averaging over a region of effective radius 
170 AU. 

For this region of the northern jet, Itoh et al. (2000) have obtained a spectrum in the 
wavelength interval 1.05 to 1.82 yum and identified eight [Fell] lines. These lines would be 
less affected by intervening extinction than the optical [S II] lines. Also, as noted by Itoh et 
al. (2000), most critical densities (ratio of radiative to collisional rates) for [Fell] are higher 
than those of [SII]. Consequently, the [Fell] lines trace potentially gas at higher densities. 

We have extended the analysis by Itoh et al. (2000) by model fitting their entire observed 
spectrum, using an [Fell] model atom having 142 energy levels with 1438 transitions and 
with atomic data from Quinet et al. (1996) and Zhang & Pradhan (1995). The absolute flux 
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calibration was kindly communicated to us by Y. Itoh. The transitions under consideration 
are identified in the energy level diagramme of Fig. 1. 

Our best fit model is shown in Fig. 2 for the spectral resolution of 280 (Itoh et al. 2000). 
These calculations assume a Full Width Half Maximum of 150 km s _1 for the [Fe II] lines (Pyo 
et al. 2002; Fridlund et al. 2004), solar chemical composition and an electron temperature 
T e = 10 4 K. Little is known about the abundances in the jet, but the assumption of a solar 
iron abundance seems not too unreasonable (Beck-Winschatz & Bohm 1994; Bohm & Matt 
2001). The absence of lines from highly excited states in the observed spectrum limits the 
temperatures to significantly below 3 x 10 4 K, a regime in which the final results are not 
critically sensitive to the assumed temperature. 

For optically thin emission (verified by the model calculations), the line ratio /i.64^m/-^i.25 ; 
constant = 0.96. For an average interstellar extinction law (Rieke Sz Lebofsky 1985), the 
observed line ratio of 1.6 implies a visual extinction, Ay, of 5.72 magnitudes. This value 
is consistent with the absence of detectable Paschen line emission, based on the Ha-flux of 
Fridlund et al. (2004), F a = 4 x 10~ 15 erg cm" 2 s -1 , and on CaseB recombination (Hum- 
mer & Storey 1987). This Ay- value is also an upper limit, since for higher values, the 
Pj3 A 1.2818 fim line would be readily seen in the spectrum of Itoh et al. (2000). Similar 
applies to P7 A 1.0938 /im. 

With the extinction fixed, the fitting of the observed spectrum results in a local jet 
density, formally, as logn e = 4.88 ±0.02, i.e. (7.6±°j) x 10 4 cm , for la observational errors 
(Itoh et al. 2000). These densities in the northern jet are much higher than what can be 
determined for the southern jet. 

The source size of this model is 60 AU (jet diameter = 0'.'4), the column density is 
7V [FeII] = 2x 10 15 cm" 2 and the total [Fell] cooling rate is L [FcU] = 7x 1O" 4 L (2.7 x 10 30 erg 

2.2.3. The radio jets 

On the O'.'l to arcsec scales, Rodriguez et al. (2003a) have recently confirmed the binarity 
of the L 1551 IRS 5 jet also at radio wavelengths (3.5 cm). The radio jets appear well aligned 
with their optical counterparts which become detectable only further 'down- stream', because 
of heavy extinction. Given that the northern jet is the dominant one at most wavelengths, 
it might seem enigmatic that the 3.5 cm emission from the southern jet is stronger by a 
factor of about two. For free-free emission, a relatively mild, by a factor of \pl, variation 
of the electron density, integrated along the line of sight (S^frce-frcc oc /(x e nn) 2 dZ), could 
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accomplish this. The southern jet does not conform with the generally adopted idea about 
the disk-jet geometry, as its direction appears to deviate substantially from the expected 
orthogonality (Rodriguez et al. 2003a). 

2.3. Properties of the binary protostar 

On the basis of high resolution VLA observations at millimeter wavelengths, Rodriguez 
et al. (1998) resolved the central source IRS 5 into two components. The marginally resolved 
emission peaks have an elliptical appearance and most likely represent circumstellar disks 
around each of the components in a protostellar binary system. The protostars themselves 
are, of course, not directly detected at the wavelength of 7 mm. 

2.3.1. The photospheric spectrum 

The spectroscopic observations of the nebulosity HH 102 (S239) by Mundt et al. (1985) 
revealed the reflected photospheric spectrum of the deeply embedded object IRS 5. The 
wavelength range was extended by Stocke et al. (1988), who concluded that the data were 
consistent with stellar spectral types of giants (luminosity class III) ranging from G2 (in 
the blue spectral range) to K0 (in the green spectral region), with an uncertainty of two 
sub-classes. In addition, the absorption line depths were indicative of a surface gravity even 
lower than that of supergiant stars (luminosity class I). Taken to the extremes, this would 
allow for a considerable range in photospheric temperatures, viz. from 4300 K (K2 1) up to 
5600 K (GO III), see, e.g., Cox (1999). In our analysis below, we will adopt a slightly less 
conservative range, viz. T e fj = (5100 ± 200) K, which covers the spectral types K0III to 
GO III and which includes the temperatures of supergiants of earlier spectral type. Stocke 
et al. (1988) interpreted this drift in spectral type, in combination with the low gravity, as 
evidence for an FU Orionis (FUOR) type of disk around IRS 5. 

As already noted by Mundt et al. (1985), and reinforced by Stocke et al. (1988), 
the Balmer line profiles exhibit P Cygni structure, indicative of a wind velocity of about 
440 km s -1 , essentially identical to the Ha line width observed by Fridlund et al. (2004) in 
emission from the northern jet. 
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2.3.2. The radiative luminosity 



For a distance of 150 pc, detailed and self-consistent fitting of the entire observed spectral 
energy distribution (SED) of IRS 5, using a two-dimensional radiative transfer model for a 
disk structure, led to the determination of the total luminosity, L tot = 40 L (White et al. 
2000). This value is larger by nearly 40% than the calorimetric luminosity, L ca i = 30 L , 
obtained from direct integration of the SED, owing to photon escape in the low-density 
polar directions. The model also correctly reproduces observed spatial intensity profiles and 
interferometric visibilities, lending further confidence in the luminosity estimate by White 
et al. (2000). 



From multi-epoch radio observations of the proper motions of the binary source IRS 5, 
Rodriguez et al. (2003a) have estimated the dynamical mass of the system as 0.1 < m*, j+ 
Si m D.i ~ 1.2 M Q . The mass of the disks X]j m D,i — 0.1 Mq, with the northern disk being 
twice as massive as the southern one (Rodriguez et al. 1998). 



For observationally derived values of 9 (45° to 60°, Sect. 2.2.1), the velocity of the 
northern jet, > 600 km s" 1 , is comparable to the escape velocities of main-sequence stars, 
given by 



with obvious notations. Since the combined mass of the binary amounts to about 
IMq, the total luminosity should not exceed 1L Q , if the stellar components were in the 
main-sequence. This is inconsistent with observation (L tot = 40 L & ) and attests to the pre- 
main-sequence nature of the objects. One notes that this luminosity, on the one hand, is 
much too low for a typical FUOR, but also much too high for a T Tauri star, on the other. 

It is illustrative to try to place IRS 5 into the H-R diagram, as shown in Fig. 7 of 
Stahler et al. (1980), where it would be somewhere near the top of the curve labelled 'gas 
photosphere'. This would also be consistent with the age of the large scale molecular outflow, 
viz. ~ 10 5 yr (Snell et al. 1980; Padman et al. 1997). In the figure by Stahler et al. (1980), 



2. 3. 3. The dynamical mass of the system 



3. 



Protostar models 




(1) 
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the evolutionary tracks for the protostellar 'dust' and 'gas photosphere', respectively, are 
separated. However, unlike their case of isotropic emission, the protostellar photosphere of 
IRS 5 can be viewed in reflection, since the dust 'shell' is not optically thick in all directions. 

At the time Stocke et al. (1988) wrote their article, the binary nature of IRS 5 was not 
established. As an alternative to their suggestion, we shall below explore the possibility that 
the scattered light observed by Mundt et al. (1985) and Stocke et al. (1988) is due to the 
combined spectrum of two protostellar photospheres. 



At present, only the total luminosity, L tot) of the IRS 5 system is determined obser- 
vationally. Assuming that the binary is protostellar in nature and, as such, derives its 
luminosity mostly from mass accretion processes, the total luminosity is the sum of the 
accretion luminosities of each of the members of the system, i.e., with common notations, 



We assume that the stellar contribution to L tot = 40 L Q (White et al. 2000) is provided 
by two objects of total mass, M tot = 1.1 M (Rodriguez et al. 2003b), with a minimum 
mass of about 0.1 M for one of the components. Thus, for the examination of Eq. 2, we 
can limit the mass range for the individual masses to Mj e [0.1, 1.0] M & . This is 'fortunate', 
since in this mass interval the mass-radius relation, yielding the ratio (M/R)i for a given 
mass accretion rate M acCt i, is insensitive to the details of the accretion processes and physical 
boundary conditions (see Stahler 1988; Palla & Stahler 1992, 1993, and Fig. 3), provided 
these objects are in their deuterium burning stage. 

The cited references provide mass-radius relations, R(M), for a few mass accretion rates. 
For other values, we obtained the mass-radius relations by interpolating in between these 
rates (see Fig. 3). 



The effective temperatures of the protostellar photospheres are obtained from (Stahler 
1988; Palla & Stahler 1993) 



3.1. 



The mass-radius relations and accretion luminosities 




(2) 



3.2. The surface luminosities and effective temperatures 
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(3) 



where the surface luminosity, L surf j, is the sum of the radiative and convective contri- 
butions, given by 



with Lq = 0.153 Lq, and where we have approximated L corw ,i with the luminosity during 
full deuterium burning, Ld,i- For accretion rates different from those given by Stahler (1988), 
L^i- values were obtained by interpolating the published Lrj(M)-curves. Since our primary 
objective is to obtain some reasonable estimates of T e g t i, we did not bother to attempt 
extrapolating the LD(M)-curves beyond those given by Stahler (1988). Instead, for 1.1 x 1CT 5 
to 1.5 x lO~ 5 M yr _1 we used the curve for 1 x 10~ 5 M Q yr _1 , a procedure which will not 
affect our conclusions below. 



We examine numerically Eq. 2 on the intervals Mi G [0.1, 1.0] M and M aC c,i G [2, 15] 10" 6 M yr" 1 , 
within which the parameters for both binary components are allowed to vary on the adopted 
grid, SM = 0.025 M and 5M acc = 1.0 x 10 _6 M o yr _1 , respectively. This discretisation 
introduces a 'fuzziness' on the boundary condition L tot = 40 L Q , which we estimate from 



(AL/L) acc ~ ±[{dM/M) 2 + {dM/Mf + {dR/R) 2 ] 1/2 as about ±15% (±5.8 L & ). This is 



comparable to the observational uncertainty. As an additional constraint, the effective tem- 
perature of the more luminous component is bounded by the interval T e ff.j G [4900, 5300] K 
(see Sect. 2.3.1). 

We find a couple of islands of formally acceptable solutions, shown in Fig. 4. There, 
the explored parameter space for the variables mass, accretion rate, accretion luminosity 
and effective temperature is depicted. Most solutions select the more massive star (the 
'primary') as also the more luminous one. This is shown by the larger area encompassed 
by the full-drawn white curve. The adjacent dashed lines outline the region, where the 
corresponding secondaries are situated. 

A smaller number of other solutions were also found, where the less massive star (the 
'secondary') is the more luminous one, because of a much larger mass accretion rate (and, 




3.3. Physical parameters of the binary components 
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hence, becomes the primary in the commonly accepted sense). These are shown by the two 
separated islands. 

Our preferred solution is shown in Fig. 4 by the two black dots, where the larger one 
signifies the primary (both more massive and more luminous). The reasons for this selection 
will become apparent below. The physical parameters for the protostellar primary of IRS 5 
are: M x = 0.8 M , M acc ,i = 6 x RT 6 M yr" 1 , T cfM = 4.9 x 10 3 K, log#i = 3.1 and 
^esc,i = 270 km s" 1 . For the secondary, the corresponding values are M 2 = 0.3 M , M acCi 2 = 
2 x lO^Moyr -1 , T effj2 = 5.8 x 10 3 K, log# 2 = 3.2 and w csc , 2 = 235 kms~\ respectively. The 
secondary contributes 25% to the total luminosity and the derived values of log g correspond 
to those of normal giants (luminosity class III, cf. Sect. 2.3.1). 

3.3.1. Photospheric emission 

The parameters derived in the previous section can be used to predict the photospheric 
spectra of IRS 5. Simple estimates of the resulting monochromatic luminosities indicate that 
the photospheric emission of the hotter secondary should be comparable in intensity to that 
of the cooler primary in the blue spectral region, but weaker already in the green (and for 
longer wavelengths). This is verified in detail when examining theoretical stellar atmosphere 
models. We use the NextGen models of Hauschildt et al. (1999) for the appropriate effective 
temperatures, the closest available values of the surface gravity (logy = 3.5) and solar 
chemical composition. 

The individual model spectra and their sum are shown in Fig. 5, whereas in Fig. 6, the 
normalized composite model spectrum is displayed together with the observations by Stocke 
et al. (1988). The latter models have been 'spun up' to the escape velocities of the individual 
components prior to adding them together. These velocities must be regarded as upper limits 
to real and observable rotation speeds (break-up and viewing geometry, respectively), but 
these values are consistent with the limited spectral resolution of the observations of Stocke 
et al. (1988). 

Rotation at, or rather close to, break-up in the stellar equatorial regions would also 
naturally explain the extreme supergiant characteristics of some of the absorption lines, 
since in this case, logy corresponds by definition to zero gravity. 

For the interpretation of the nature of IRS 5, the implications of the FUOR-disk and 
protostellar photosphere scenarios are very different. High resolution spectroscopy in the 
optical should enable us to distinguish between these alternate models and, furthermore, 
potentially provide the opportunity to study stellar surfaces during an evolutionary phase 
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which has not previously been accessible to direct observation. 



3.4. Outflow models 

To put the observed and derived properties of the binary protostar and the binary jet 
into context, we will make use of theoretical models of outflows. Specifically, the theory of 
magnetocentrifugally driven flows (x-winds) has been worked out in considerable detail and 
presented by F. Shu and co-workers in a series of papers. 



3.4-1- x-wind velocities 

In terms of the mass-radius relation (Sect. 3.1), the expression for the terminal wind 
velocity by Shu et al. (1994, their Eq. 4.13a) can be recast into 



(5) 



where r x is a factor, such that r x times the stellar radius R is the radial distance of the 
location of the x-point, assumed to be close to the radius of corotation of the stellar surface 
and the circumstellar Keplerian disk, so that r x is of order unity. The star is assumed to 
rotate near break-up. 

As explicitly indicated, Eq. 5 represents an upper limit to the wind velocity, with equality 
reached at infinity (Shu et al. 1995). The angular momentum parameter, J, takes values 
> 1.5 and is likely not to exceed ~ 10. For the mass-radius relation corresponding to the 
accretion rate of 6 x 10~ 6 M® yr _1 (Fig. 3), graphs of Eq. 5, with r x = 1.0, are shown in Fig. 7. 



3-4-2. x-wind f -factors and densities 

For a single stellar mass and a particular choice of parameters (M = 0.5 M & , R = 4 R Q , 
r x = 1.0, M w = lO~ 6 M yr _1 ), Najita & Shu (1994) presented detailed numerical results, 
covering the range in J from 2.0 to 7.8, with accompanying magnetic field strengths from 
1.5 to 8.3 kG. Corresponding /-factors vary from 0.4 to 0.1, where / is defined as the ratio 
of the mass loss to the mass accretion rate. 

The collimation of, initially wide-angled, x-winds into narrow jets has been addressed 
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by Shu et al. (1995). These authors also make a theoretical prediction of the density profile 
across the jet, i.e. in terms of the jet radius, the density scales approximately as p(r) oc r~ 2 . 

4. Discussion 

4.1. The deuterium abundance 

The results presented for the protostar binary IRS 5 are based on the mass-radius re- 
lations during deuterium burning as provided by the models of Stahler (1988) and Palla & 
Stahler (1992). These were calculated for the deuterium abundance relative to hydrogen, 
[D/H] = 2.5 x 10" 5 , and such, sensitive to the assumed [D/H]. 

Recent observations indicated the significantly lower value of about 1.5 x 10~ 5 in the 
solar neighbourhood and, furthermore, with a considerable dispersion (Sonneborn et al. 2002; 
Moos et al. 2002; Steigman 2003). However, Linsky & Wood (2004) recently announced that 
they could reconcile the observational data with a value of [D/H] = (2.3 ±0.4) x 10~ 5 within 
lkpc of the Sun, in which case the results of the calculations by Stahler (1988) and Palla & 
Stahler (1992), with regard to the assumed deuterium abundance, remain valid. 

4.2. Single-star theory and binary system 

We obtain a binary mass ratio of q = MijM\ ~ 0.4, which is close to the peak of 
the observed distribution for solar-type stars in the solar neighbourhood (Duquennoy & 
Mayor 1991). These authors did also not find any statistically significant dependence of the 
g-distribution on binary period (hence, binary separation). 

Both the photospheric spectrum and the jet momenta are not easily reconcilable with 
an equal-mass binary, accreting at the same rate. Also the higher disk mass of the primary 
would heuristically be compatible with the higher M acCi i (M acc oc S a _disk, e.g. Frank et al. 
1985), since both disks have the same size (about 10 AU, Rodriguez et al. 1998). 

Above, we have tacitly assumed that the components of the protostellar binary can be 
modelled with the theory of single objects. The relatively large separation (43 AU/ cos binary > 
a crit ; Rodriguez et al. 2003a; Hale 1994) justifies this assumption. Here, a crit ~ 30 — 40 AU 
is the 'critical separation', beyond which observed solar-type binaries cease to be coplanar 
(Hale 1994), and hence can be assumed to have developed independently. 
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4.3. The binary jet-driving IRS 5 and x-wind models 

Each component of the protostellar binary drives its own jet. The (total) wind mass 
loss rates of Sect. 2.1 are therefore likely upper limits to the individual loss rates. Comparing 
with the result of Sect. 3.3, the ratio of the mass loss to accretion rate for the primary is 
estimated at fi = (M w /Macc)i < 0.13 (< 0.33), where the value in parentheses refers to the 
HI rate. Correspondingly, for the secondary, f% < 0.4 (< 1.0). 

When compared to x-wind models, these values are consistent with high J (~ 8 — 10) and 
low J (~ 2), respectively (Najita & Shu 1994, see their Table 5). The inclination corrected 
jet velocities conform with this scenario (see Fig. 7). Based on our analysis, we identify the 
driver of the fast and well-aligned northern jet with the protostellar primary (IRS 5-N) and, 
consequently, the source of the much slower southern jet with the secondary (IRS5-S). 

Unless the disks are occasionally rejuvenated by the significantly more massive cir- 
cumbinary reservoir (Osorio et al. 2003; Fridlund et al. 2002, and references therein), the 
current level of the mass loss through the wind could be maintained for at most another 
M D ,i/Macc,i ~ 0.06 M /6 x lO^Moyr^ 1 = 10 4 yr. This would be only about one tenth of 
the age of the large scale CO outflow, but comparable with the time scale (~10 3 yr) of the 
molecular outflow close to the source (Fridlund & Knee 1993). Given the large difference 
in their momentum rates, these two flow phenomena appear only indirectly related. The 
scenario which emerges (Fridlund et al. 2004) is that the molecular flow lifts off the disk at 
some AU distance (disk-wind, e.g., Konigl & Pudritz 2000), whereas the atomic flows/jets 
originate much closer to the protostellar surface (x-wind). 

The results of Sect. 2.2.2 can be used to examine the density behaviour across the 
northern jet. The observed difference in density on different spatial scales hints at a non- 
flat density distribution in the jet. These estimates of the density are based on optically 
thin line emission and, as such, represent averages over volume, i.e. n = Jn dV/ f dV. We 
assume that the density distribution perpendicular to the jet axis can be approximated by 
a power law, n(r) oc r a . If we examine a sufficiently small section of the cylindrical jet far 
away from the source, so that the density distribution perpendicular to the jet-axis does not 
depend on the height-coordinate (^-coordinate) of the cylinder, the integrals become trivial 
and n = 2n(r ){R/r ) a [l - (r /R) a+2 ] /{a + 2), a ^ -2 , and n = 2n{r ){R/r ) a In (R/r ), 
a = —2 , where ro -C R is some fiducial radius close to the central jet axis. This recovers 
the well known result that the average density is dominated by the largest scales. In either 
of these two cases, the power law exponent is given approximately by a ~ A log n/ A log/?, 
as usual. Our line data represent post-shock values for gas which has been compressed by 
an, in general, unknown amount. However, for a planar cross section, the compression is 
constant over the shock surface and from the [Fe II] and [S II] data we then obtain formally 
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ol = —2.3^2 for the observed, and resolved (Fridlund et al. 2004), scales -Rpeii] = 0"2 and 
R\su] = 1"0. This observationally determined value of the power law exponent is in reasonable 
agreement with the theoretical one for x-winds, viz. a — > —2 as z — > oo (Shu et al. 1995). 

This circumstance motivates the direct comparison with the predictions by the x-wind 
model. With the formalism of Shu et al. (1995) and for the parameters of IRS5-N, we 
estimate that the central density (at r = R x ) is p x < 1.8 X 10 -12 (4.6 x 10 -12 ) gem -3 , with 
the same convention regarding the mass loss rate as before. At the deprojected distance from 
the protostar of about 400 — 600 AU and at the radial distance from the jet axis r = 30 AU, 
we estimate from Fig. 3 of Shu et al. (1995) that the jet density should be of the order of 
6 x 10 -20 (15 x 10 -20 ) gem -3 , which translates into a neutral particle density of 4 x 10 4 
(9 x 10 4 ) cm -3 of unshocked gas. As already remarked above, the observationally derived 
value (~ 7.5 x 10 4 cm -3 ) refers to the post-shock electron density, assuming an ionisation 
fraction of unity. The comparison of these two values suggests that the compression of the 
shocked gas is just offset by the fractional ionisation of the neutral jet flow. This seems 
reasonable, as the former is, on the average, a factor of about 80 (Hollenbach et al. 1989), 
whereas the latter is a few percent (Shang et al. 2004). 

Shang et al. (2004) have recently applied the x-wind model to the radio jet emission 
from IRS 5. These authors focus on the southern jet because of its apparent higher radio flux 
density. However, this circumstance is not reflected at other wavelengths, including X-rays, 
optical, infrared, millimeter and the radio regime (for the latter, see Rodriguez et al. 1998), at 
which the northern source and jet dominate. In addition, the flux difference is small enough 
that it could be easily absorbed in the model by Shang et al. (2004) and we conclude that the 
x-wind theory provides a viable model which is capable of explaining various independent 
pieces of observational evidence. A final and more decisive test of the theory would include 
the direct measurement of the magnetic field and the protostellar rotation rate. 



5. Conclusions 

Our main conclusions from this work can briefly be summarised as follows: 

• Interpreting recent observational achievements within the framework of the theory of 
protostellar structure and evolution allows the derivation of the physical properties of 
the individual components of the binary protostar L 1551 IRS 5. 

• We offer an alternate scenario to the FUOR-disk hypothesis to explain optical spectro- 
scopic data of the scattered light from IRS 5, which can be put to test with currently 
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available observing facilities. This offers potentially the unique opportunity to directly 
observe the surfaces of accreting protostars. 

• The derived properties of the binary, in combination with their observed mass loss rate, 
lead to results which are consistent with detailed predictions by the theory of x- winds. 



We thank the anonymous referee for her/his thoughtful comments. RL enjoyed inter- 
esting discussions with Mary Barsony, Goran Olofsson and Frank Shu. We are grateful 
to Francesco Palla for making available to us the mass-radius relations for different shock 
boundary conditions. 
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Fig. 1. — The thirteen lowest energy levels of [Fell] are displayed. The 19 transitions marked 
in the spectrum of Fig. 2 are shown by the lines connecting various levels. The two strongest 
lines in that spectrum are shown by the thick lines. 
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Fig. 2. — The observed and synthetic spectra of the northern jet in the photometric J- 
and H-bands. The black curve corresponds to the flux densities calculated from the line 
flux observations of Itoh et al. (2000, Table 1, ID=Eb, R x = 280), who identified 8 [Fell] 
lines, the transitions of which are labelled in green (cf. Fig. 1). Other line fluxes were 
estimated from their Fig. 5. The fitted model is shown in red: whereas the matching in 
the J-band is excellent, some of the weaker lines in the H-band appear too strong in the 
model. This 'mis-fit' could be due to insufficient background correction. The OH lines are 
particularly intense in this wavelength region, as is illustrated by the blue spectrum, showing 
the arbitrarily scaled telluric OH emission (Rousselot et al. 2000). Since X-ray emission has 
been detected from this jet position (Favata et al. 2002; Bally et al. 2003; Fridlund et al. 
2004), we tentatively identify the emission feature at 1.08 /xm as the recombination line 
HeIAl.0830 ( 3 P° 2p - 3 S2s). 
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Fig. 3. — Mass-radius relations of protostars during deuterium burning. The curves drawn 
with thick lines are adopted from the literature, with our interpolations shown by the thin 
lines. The parameters next to the curves are the mass accretion rates, with the data for 2 x 
10" 6 , 5 x 1(T 6 and 1 x 10" 5 M yr" 1 taken from Stahler (1988) and those for 3 x 1(T 5 M yr" 1 
from Palla & Stahler (1992). 
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Fig. 4. — The explored parameter space in the variables mass accretion rate, accretion 
luminosity, effective temperature and protostellar mass. Here, the radii of the protostar 
models are implicit and presented separately in Fig. 3. The bar to the right of the graph 
provides the colour coding for the effective temperature. Areas, thought best to represent 
the protobinary L 1551 IRS 5 are encompassed by the white lines, where the full-drawn lines 
indicate the positions of the primaries and the dashed lines those of the companions. The 
smaller and clearly separated areas correspond to pairs in which the less massive star is the 
more luminous object. The black lines are meant to facilitate the viewer's orientation along 
the axes. Our best solution for the IRS 5 binary is shown by the fat black dots, with the 
bigger one identifying the primary (IRS5-N, see the text). 
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Fig. 5. — Stellar model atmospheres from Hauschildt et al. (1999) for log g = 3.5, solar com- 
position and the indicated T e g, representing (approximately) the preferred binary solution 
shown in Fig. 4. The hotter secondary is comparable in brightness to the cooler primary in 
the blue spectral regime, resulting in the relatively shallower G-band near 4300 A than what 
would be expected on the basis of the green spectrum, where the primary is dominating the 
photospheric emission. The upper curve is the sum of the spectra in the two lower panels. 
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Fig. 6. — The normalised observed optical spectrum of IRS 5 (Stocke et al. 1988) shown 
by the thick lines and the composite binary spectrum (thin lines). The latter is based on 
the models shown in Fig. 5, but with their spectra rotationally broadened according to their 
escape velocities, prior to summation (see the text). Several features, such as H54101 and 
H/5 4861 in particular, seem not to match very well, which can be ascribed to contaminating 
emission from the observed nebula HH 102 (Stocke et al. 1988). 
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Fig. 7. — Various velocities versus protostellar mass. The deprojected radial velocities of the 
jets from L 1551 IRS 5 are shown by the shaded regions labelled Jet-N and Jet-S, respectively. 
The dots denote the positions of the corresponding driving sources, viz. the primary IRS 5-N 
and its less massive companion IRS5-S. These are shown on curves of the escape velocity, 
corresponding to their mass accretion rates of 6 x 10~ 6 and 2 x lO~ 6 M yr _1 , respectively. 
For comparison, the escape velocities for the main-sequence (ms) are also shown by the 
dotted line. The curves with the upper limit symbols are the loci of the asymptotic x-wind 
velocity for different values of the angular momentum parameter J and the accretion rate of 
6 x 10~ 6 M Q yr" 1 . To reach velocities as high as those of Jet-N, large values of J are required, 
implying /-values of 0.3 or smaller (see the text). 



